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We study the situation where two point like mirrors are placed in the vacuum state of a scalar 
field in a two-dimensional spacetime. Describing the scattering upon the mirrors by transmittivity 
and reflectivity functions obeying unitarity, causality and high frequency transparency conditions, 
we compute the fluctuations of the Casimir forces exerted upon the two motionless mirrors. We use 
the linear response theory to derive the motional forces exerted upon one mirror when it moves or 
when the other one moves. We show that these forces may be resonantly enhanced at the frequencies 
corresponding to the cavity modes. We interpret them as the mechanical consequence of generation 
of squeezed fields. 

PACS: 03.65 - 42.50 - 12.20 



INTRODUCTION 



' The vacuum fluctuations of the electromagnetic field manifest themselves through macroscopic Casimir forces [[l] [| . 
7-H , These forces result from the radiation pressure exerted upon boundaries by the scattered fluctuations which depends 
upon the instantaneous and local value of the stress tensor. Consequently, they are fluctuating quantities §. More 
precisely, the forces exerted upon the two mirrors of a Fabry-Perot cavity have to be considered as random variables. In 
the present paper, we compute the associated correlation functions. For the sake of simplicity, we study the situation 
where two point like mirrors are placed in the vacuum state of a scalar field in a two-dimensional (2D) spacetime. 

As illustrated by the Langevin theory of Brownian motion ||, any fluctuating force has a long term cumulative 
effect. This cumulative force can be derived from linear response theory 0. The fluctuations of Casimir forces 
thus imply that mirrors moving in the vacuum must experience systematic forces. In a previous paper J?j], we have 
computed such a motional force for a single mirror in the vacuum state (of a scalar field in a 2D spacetime). At 
the limit of perfect reflectivity, we found a dissipative force proportional to the third time derivative of the mirror's 
i position q 



which corresponds to a linear susceptibility at the frequency to 

xH = pi (ib) 

D7TCT 

(from now on, we use natural units where c = 1; however, we keep h as a scale for vacuum fluctuations). This 
damping force for a single moving mirror is connected to the Casimir force (mean force between two mirrors at rest) 
since both result from a modification of the vacuum stress tensor ||. Actually, expression (la) identifies with the 
linear approximation of the force computed using the techniques of quantum field theory As required by 

Lorentz invariance of the vacuum state ]l2]| , the damping force vanishes for a motion with a uniform velocity. In the 
case of a uniform acceleration, the mirror is submitted to the same fluctuating field as if it were at rest in a thermal 
field (ll| so that it experiences also a zero motional force. 

In the present paper, we compute the explicit expressions of the forces exerted upon one mirror due either to its 
own motion (in presence of the second mirror) or to the motion of the other. These expressions, obtained from linear 
response theory, are valid in a first order expansion in the mirrors' displacement without restriction on the motion's 
frequency. 

A problem in any calculation of vacuum induced effects is to dispose of the divergences associated with the infinite- 
ness of the total vacuum energy. This problem can be solved by assuming that the boundaries are transparent at high 
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frequencies. Using a scattering approach where the mirrors are described by transmittivity and reflectivity functions 
obeying unitarity, causality and high frequency transparency conditions, one obtains a regular expression for the mean 
force between two mirrors [[|. The same approach also provides directly a causal motional force in the single mirror 
problem whereas the non causal expression (1) is recovered as an asymptotic limit for a perfectly reflecting mirror [M. 

In the present paper, we use this approach to study the motional forces in the two mirrors problem. First, we 
compute the correlation functions characterizing the fluctuating Casimir forces exerted upon the two motionless 
mirrors. Then, we use the linear response theory to derive the susceptibility functions associated with the motional 
forces. In order to obtain these functions, we use some analytic properties of the correlation functions which are 
analysed in Appendix A. We check that our results are consistent with already known limiting cases: the static 
Casimir force J3| (limit of a null frequency), the one mirror problem (limit where one mirror is transparent) and 
the limit of perfect reflection [ |To| |. 

The expressions obtained for perfectly reflecting mirrors correspond to a damping force analogous to equation (1) 
with two differences. First, the response is delayed because of the time of flight from one mirror to the other. Second, 
the motional modification of the vacuum fields is reflected back by the mirrors. The resulting interference between 
the different numbers of cavity roundtrips gives rise to a divergence of the susceptibility functions. For partially 
transmitting mirrors, these functions are regular and describe a resonant enhancement of the motional Casimir force, 



which may be large when the Fabry-Perot cavity has a high finesse 14 . A resonance approximation is used to evaluate 
the susceptibility functions in this case. 

It is known that mirrors moving with a non uniform velocity squeeze the vacuum fields |l5| , |l6| . This requires that 
energy and impulsion be exchanged between the mirrors and the fields. The motional Casimir forces thus appear 
as a mechanical consequence of this squeezing effect. In Appendix B, we discuss the motional modifications of the 
field scattering by the mirrors and we write an effective Hamiltonian describing the squeezing effect as well as the 
mechanical forces upon the mirrors. 

NOTATIONS 

Any function f(t) defined in the time domain and its Fourier transform f[uj] are supposed to be related through^] 

/(*) = / ^fW™ (2a) 
In the following, some functions of time will be expressed as integrals over two frequencies, with the notation 

m = \%\%^ t - t ^'\ (2b) 

Comparing with (2a), one gets the equivalent expression 

f[u)= j^f[u-u',J] (2c) 

In a 2D spacetime (time coordinate t, space coordinate x), a free scalar field is the sum of two counterpropagating 
components <p(t — x) + ij)(t + x) which will be written as the two components of a column matrix 

$ (t) - ( ^ ~ X) 

The Fourier transform of the column <f> x {t) is related to the standard annihilation and creation operators 

corresponding to the two propagation directions 

*.m -($£:) -^m 



1 The notation used in the original paper for Fourier transforms has been changed to a more convenient one. 
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where the abbreviated notation $ stands for the values of $> x evaluated at x = and 

V = 



1 

-1 



The energy and impulsion densities correspond to two counterpropagating energy fluxes 
e x {t) = <p'(t ~ x) 2 + ip'(t + x) 2 p x {t) = ip'(t - xf - i//(t + xf 
They may be written as integrals over two frequencies (see eqs 2) with 



e x [u>, to'} = iojioj' Tr [^[cjJ"^^ 



e x \uj,uj 



Tr stands for the trace operation and X T for the transposed of X. So, their mean values in a quantum state are 
functions of the covariance matrix 

{^ x [lu]^ xI [lo'] t ) = e^ x ($[w]$[w'] T )e i '*'' V 

For a stationary state such as the vacuum, the covariance matrix depends only upon one frequency parameter 

($M*[u/] T ) = 2m$(w + w')c[uj] (3a) 

It will be useful to write it in terms of the anticommutators which characterize the field states and of the commutators 
which do not depend upon the states 

c[uj] = c+[lu] + c_[w] (3b) 



c+ M = 



c[u>] + c[—ui] T 



c\u>\ — c — u>\ 



in 

4w 



/ is the unit matrix. The vacuum state corresponds to 



c[w]=J0(w) — 



c + [w]=Ie(w)— e(uj) = 9(uj) - 



(4) 



SCATTERING UPON MOTIONLESS MIRRORS 



The scattering of the field upon a partially transmitting mirror at rest at q is described by 

The S— matrix is supposed to be real in the temporal domain, causal and unitary 

S[-u] = S[uj}* 

S[uj] is analytic (and regular) for Imcj > 



S[u]S[l 



It - 



1 



(5) 

(6a) 
(6b) 
(6c) 

(6d) 

This assumption will allow us to regularize the ultraviolet divergences associated with the infinitencss of the vacuum 
energy ||. A mirror perfectly reflecting (s = and r = —1) at all frequencies does not obey this condition and it 
will be better to consider the perfect mirror as the limit of a model obeying the transparency condition (for example 
a mirror perfectly reflecting at frequencies below a reflection cutoff) . 

We study now the situation where two mirrors are placed at rest in the vacuum at positions q\ and q 2 (see Figure 
!)• 



Finally, the mirror is supposed to be transparent at high frequencies 

S\lu] — > / for lo — > oo 



Pi, 



<p 



out 



FIG. 1. Two mirrors scatter the two counterpropagating fields. The subscripts 'in', 'out' and 'cav' refer respectively to the 
input, output and intracavity fields. 
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As the final results will depend only on the distance q — g 2 — qi between the two mirrors, we consider from now on 

92 = g qi = — 2 (7) 

The scattering of the field by the mirrors (see Figure 1) is described by two scattering matrices given by equations 
(5) 

VcavM A _ / Si[(J\ n[uj]e luJq \ ( <Pin{u] 

^outH J \ ri [u)]e-^i siM J \i> cav [u] 

VoutM \ = / S2H T 2 [uj\e~ lUq \ ( <^ C av[H 

V'cavH J \r 2 [uj}e 1 ^ s 2 [uj} J \ Tp in {u] 

As usual in the computation of a Fabry Perot cavity, one can solve these equations to express the outcoming and the 
intracavity fields in terms of the input ones 

$outM = SM^inM <$> cav [u] = R[u]<f> in [u] 

The global scattering matrix S[lu] and the resonance matrix R[uj] are given by (all the reflectivities and transmittivities 
are evaluated at frequency lo) 

q , ,_ 1 ( s lS2 r 2 e-^i + tferie^ \ m 

l J d[u] \rie-*"* + dir 2 e 1 ^ sis 2 J w 



1_ ( s i 

d[w] V si^e 1 " 9 s 2 

di[uj] = si[w] 2 - ri[uj] 2 d 2 H = S2M 2 ~~ r 2M 2 

= 1 - r[uj]e 2lujq r[w] = ri[w]r 2 [o;] 



One checks from conditions (6) that the matrices 5 and i? are real in the time domain, that they are retarded response 
functions, that the matrix S is unitary and that the Fabry-Perot is transparent at the high frequency limit 

S[-u] = S[w}* R[-u] = R[(j]* (10a) 

S[u>] and R[ui] are analytic (and regular) for Imw > (10b) 

S[v]S[u]i = 1 (10c) 

S[u] -> I R[uj] — > I for w -> 00 (lOd) 



CASIMIR FORCE BETWEEN TWO MOTIONLESS MIRRORS 

As usual in a local formulation of the Casimir effect jl7| , the forces F\ (t) and F% (t) acting upon the two mirrors 
can be deduced from the stress tensor evaluated at both sides of the mirrors 

Fi(t) = if[ n (t - qi ) 2 - ip' c&v (t - qi ) 2 + ^ ut (i + gi) 2 - < av (t + 9 i ) 2 

The forces are quadratic forms of the vacuum fields as the energy or impulsion densities, so that they can be expressed 
as integrals (2) over two frequencies with (using equation 7) 

F x [uj,u}'] = exp (-^(v + u')q\ iuid (^[^^[a/] - <A»vM9W [<«>']) 

+ exp f^(uj + uj')qj iujitjj' (V'outMV'out [<*>'] - VwM^cav [<*>'] ) 

F 2 [u!,Lli'} = exp [-^{v + UJ ')<l) iujiuJ ' (Vw[<i#cav [<•>'] - V'in M^in K] ) 

+ exp ( -(w + w')?) Ww'^cavM^cavIw'] ~ font M<Amt [w'] ) 
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Noting that the matrices 

are the projectors onto the counterpropagating components (p and tp and using the expressions of the output and 
cavity fields in terms of the input ones, one obtains 

Fi[w,u'] = iwiw' Tr [Fi[u, w']$ x [w]® x [w'] T ] (11a) 

with 

Fi[uj,uj'] = e.exp f _l( w + w ')gj (p e( - R[J] T P e .R[J\) (lib) 
+e,exp \U u + u/)q \ {S[uj'] t P^ £i S[uj] - R[u'] T P- Ez R[uj}) 

E\ = 1 £2 = -1 

The two matrices Ti obey the following properties 

JF;[u/,o/] T =Ti\JM F[oj,u']i =F[-u)',-u)] (12) 
When evaluating the mean forces in the vacuum, one obtains from equations (3,4) 

{Fi[w,w']) = 6(ui) — 2nS(u) + u') Ti F^w , -ui\ 
Using the unitarity of the S— matrix and the following expression of RR^ 



R[lo]R[uo]^ = I + Q[u] + Q[u^ (13a 
1 f rir 2 e 2iuq rie luq 



Q[ UJ } = — ( ILl^q ^t&uq ) ( 13b ) 



one shows that 

Tr^w, -u] = £l Tr (/ - R[u]R[up) = -e 4 Tr (Q[w] + Q[w]t) 
This leads to the known expression for the Casimir force between two partially transmitting mirrors 



00 



FORCE FLUCTUATIONS 



We now compute the correlation functions characterizing the fluctuations of the forces F% and F% acting upon the 
two mirrors 

Cy(t) = (^(^(0))-(^)(^) 

Inserting the expressions (11) of the forces in terms of the field operators, it appears that the correlation functions 
Cij are related to fourth order moments or second order moments of annihilation or creation operators. With the 
help of Wick's rules, the fourth order moments may be deduced from the second order ones. In simple words, the 
vacuum fields may be considered as stationary gaussian random variables and higher order statistical quantities can 
be deduced from the covariance matrix. 

Using this method, simple algebraic manipulations lead to an expression of the functions Cij as integrals over two 
frequencies (see equations 2) with 

Cij[uj,uj'] = 2wV 2 Tr [j r i [u,u']c in [u!]J r :j [uj,uj'^Ci n [u'] T ] (15) 
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This becomes for the vacuum state (see equation 4) 

h 2 

Cij[u),u)] = —9{uj)9{uj')ujLj'^ij[ijj,uj l ] 

7y[w,w'] = Tr [Filw^'^jlw,^] (16) 
In other words, the noise spectra associated with the fluctuations of the Casimir forces may be written 

dj [u] = y0(w) J" ^w'(w - w')7ij [w', w - (17) 

From properties (12), the coefficients 7y obey 

= jij[u>',oj] = = jji[-u;,-uj'] 

Equations (16-17) give the correlation functions which characterize the fluctuations of the Casimir forces exerted 
upon the two mirrors. The explicit evaluation of the coefficients 7^ in terms of the scattering coefficients is given in 
the Appendix A (see equations 27 and 21). 



MOTIONAL CASIMIR FORCES 



We study now the motional Casimir forces. These forces could be obtained by analysing the modification of the 
stress tensor associated with the mirror's motion ||. At first order in the mirrors' displacements, they can also be 
derived from the fluctuations computed for motionless mirrors by using the linear response theory We have used 
both methods and checked that they provide the same results, as it is the case in the single mirror problem fit]. Here, 
we present the linear response technique. The main steps of the other method are given in Appendix B. 

The classical motion of the mirrors corresponds to an effective perturbation of the Hamiltonian 

SH(t) = -Y i Fj(t)6q j (t) (18) 

3 

where Fj are the force operators exerted upon the two mirrors. The linear response theory provides the mean motional 
forces associated with this perturbation in terms of susceptibility functions Xij 

(8Fi(t)) = ]T fdr Xij(r)5 qj (t - r) (19a) 
3 

The susceptibility functions Xij may be deduced from the correlation functions since they are the retarded parts of 
the mean values of the force commutators £jj 

X<i(*) = (19b) 
_ [F f (0,^(0)] _ Cl3 (t)-C ]t {-t) 

The motional forces will be conveniently characterized by the susceptibility functions written in the frequency domain 
(compare with equations 1) 

(SFi[u}) =X>«M<%M (19d) 

3 

In a first step, we compute the force commutators from the correlation functions 

, Cij[u,uj'}-Cji[-uj',-uj] 
- 

Using the properties (12) and writing the covariances (15) in terms of anticommutators and commutators (see equations 
3), one shows that 
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In the vacuum state, this becomes (see equations 4 and 16) 
Writing the Fourier transform of the force commutator as 





one obtains the simple following relation with the noise spectrum (17) 

Cy[w] = 2h6(u)£ij[w] 

The noise spectrum C'ij contains only positive frequency components, as expected for the zero temperature state. 
This corresponds to the fact that the vacuum fluctuations can damp the mirrors' motion but cannot excite it fl8|| . 

In a second step, we deduce the susceptibility functions as the retarded part of the force commutators. This 
derivation relies upon the analytic properties of the correlation functions and requires a detailed inspection of the 
expression of the coefficients jij[uj,ui'] in terms of the scattering coefficients. This analysis, presented in Appendix A, 
shows that is a sum of terms which are either retarded or advanced functions of w. The susceptibility functions 

XyM are obtained by retaining only the retarded terms (see equations 19). One gets them as integrals over two 
frequencies with 

XyM = ^ (e(uh*M + e (a/)7»]) (20) 

where the coefficients 7^ are the retarded parts of the coefficients 7^ considered as an analytic function of its second 
frequency parameter 

Rr n „aikw'l . n . ,e 2iujq . n e 2iu ' q 

<y£ - +ai [u,u/]0i Kw' raH-rpr^ w' — — 

a[a;Ja[a/J d[uj\ a[uj'\ 

p2iuq p2itu'q 

+r 1 [iv']di[uj]r2[uj] , -, + ri[w]cii[a/]r2[a/] 



d[u] 1 d[w'] 

+(fl[wr+riM )( Mw re-^ + ^^)^ 

+2 -5R"dRl (21a) 

p iwq piu'q p-iuq piul' q 

7 2 >V] = -a^^a^u^'] — — - (n[a/] + n[ W ]*) (r 2 [o/] + r 2 [u]*) ( 21b ) 

ai[w,a;'] = 1 - Si[w]sj[w'] + r^w]?';^'] (21c) 
ft[w,u1 = 1 - SiM«i[w1 - ri[w]riM (21d) 

7^2 and 7^ are obtained by exchanging the roles of the two mirrors. It has to be noted that the coefficients 7^ are 
not symmetrical in the exchange of to and to' so that the function Xij M is n °t an integral restricted to the interval 
[0, u>], in contrast to the functions CV, ;[u>] (see equation 17) or 

Equations (20-21) constitute the main results of this paper. They provide the expressions of the motional Casimir 
forces (see equations 19 and 2) for two partially reflecting mirrors in a linear approximation in the mirrors' displace- 
ments. 

CONSISTENCY WITH KNOWN RESULTS 

Here, we check that the susceptivity functions (20) are consistent with already known results. This can be done in 
three limiting cases. 
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First, considering that the second mirror is transparent at all frequencies {r% = 0; s 2 — 1), one obtains 



iTiujij' 



Xu[w,a/] = — ^— (E(w)+£(u'))ai[w,w1 (22a) 
and one recovers the known susceptibility for a single partially transmitting mirror moving in the vacuum 

XiM=ihJ —u'{u-w')a x [w',u-u'] (22b) 

Then, the quasistatic susceptibilities can be computed from equation (20) 

Xij[0] = J ^Xy[w,-w] 

Xij[u>, -w] = (7^[-^, - 7y[a;, -w]) 

Using equations (21), one checks that they are consistent with the mean Casimir force (14) between two motionless 
mirrors 

Finally, one can consider the limiting case of perfectly reflecting mirrors (ri = r 2 = —1; s\ = s 2 = 0) where the 
expressions (21) may be simplified to 

Simple calculations then lead to the following expression of the motional force in the time domain 

(6F 1 (t)) = A (Sq'l'(t) + 8tf'(t - 2q) + Sq'^t - 4q) + . . . - 6q' 2 "(t - q) - 6q 2 "(t - 3q) - . . .) 

+ ^ Q^iW + H(t - 2?) + Sq'i(t - 4g) + . . . - *^(t - g) - 6q' 2 (t - 3q) - . . .) (23) 

It can be checked that this is exactly the linear approximation [l9[ | of the expression obtained for perfectly reflecting 
mirrors by Fulling and Davies [fl0|| . The terms proportional to third time derivatives appear as generalizing the 
damping force (1) for a single perfectly reflecting mirror. Now, the force exerted upon one mirror depends not only 
on its own motion but also on the motion of the other one. The response is delayed due to the time of flight between 
the two mirrors: the motional modification of the stress tensor has to propagate from one mirror to the other in order 
to exert a force on it. Moreover, the modified stress tensor is reflected back by the two mirrors. The other terms, 
proportional to velocities, are not present in the one mirror problem because of the Lorentz invariance of the vacuum 
(see the discussion in the introduction). They are associated with the existence of a static Casimir force in the two 
mirrors problem. 



RESONANT ENHANCEMENT OF THE MOTIONAL CASIMIR FORCE 

The expression (23) leads to a divergence of the susceptibility functions when evaluated at the frequencies ^ 
with k integer. Indeed, the contributions corresponding to different numbers of roundtrips give rise to a constructive 
interference at those frequencies. In other words, the Fabry-Perot can be considered as a mechanical resonator with 
resonance frequencies corresponding to the optical resonance frequencies. 

For partially transmitting mirrors, the divergences due to perfect reflection will be regularized. In the last section 
of this paper, we evaluate the resonance enhancement at the limiting case of a large quality factor. 

In this case, the reflection delays are much shorter than a roundtrip time and the reflectivity functions are smoother 
in the frequency domain than the Airy function describing the cavity modes. Replacing in (21) the denominators by 
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geometrical series, we obtain each quantity 7^ as a sum of components evolving at different frequencies: this expansion 
is analogous to the one used in the computation of the static Casimir force at the large distance approximation 

7£[w,u/] = 2ax[w,u)'] r[uj] l r[u;'] m e 2ilwq+2im "' q 

l,m>0 

+a 1 [LO,uo']f3 1 [LO,uo']r 2 [Lo]r 2 [Lo'] £ r [w] , -M^] m ~ 1 e Mw * +2W * 

l,m>l 

+ (r 1 [-u]+r 1 [u r \)r 3 [-u] ^ r[-u] l - 1 T[u/] m ■ e - 2il "«+ 2im " '? 

i>l,m>0 

+ {r 1 [-uo]+r 1 [u'])r 2 [u J '] ^ r[-uj\ l r[J\ m - 1 e~' 2a ^' llmu '« 

/>0,m>l 

+ri K]di Hr 2 [w] 2 r[w] , - 1 c MIu « 

m>l 

Z>1 m>l 



K[Lu,Lu'] = - ai [Lu,u']a 2 [Lu,u'] £ r [w] V[c/]"V( 2m 



)ojq-\-i(2m-\-l)uj' q 



721 

;,m>o 

- (nM + ri[-w]) (r 2 [w'] + ra[-w]) ^ r[-w]VK] m e- l(2i+1) ^ +l(2m+1)w ' 9 

Z,m>0 

We will denote /ijf'*^ the coefficients in this expansion 



The coefficients ^f' M ^ depend only upon the reflectivity functions. Negative values of M do not appear in the sum 
because fij^' M \u),u)'] is a retarded function of u/. Only even values of L and M appear for i = j and odd ones for 

Then, the susceptibility functions Xij M w iU be obtained through an integration (see equations 2) 

X«M = $>if ,J, °M (24a) 
do/ 

2T 



X? M) H= ^"'Ka;^ (24b) 



X g" M > [ Wj J] = ^lolo' (e( U )nlf ,M) k 0/]^+^ + £ (u/)/4j" M) [w', w ] e «"'«+ iM ^) (24c) 

It follows that the integrals Xif' with L 7^ M will contain exponentials with a rapidly varying phase and can 
be considered as non resonant terms. At the limit of perfect reflection, they provide the terms proportional to the 
velocities in equations (23) while the resonant terms Xif' with L = M provide the terms proportional to the third 
time derivatives. 

In order to obtain the behaviour of the susceptibility functions near a resonance at — with k a large integer, we 
will retain only the terms L = M (the first time derivatives have a small contribution to equation 23 compared to the 
third time derivatives in this case). In this resonance approximation, we obtain 



l>0 

,21.21) r / 

Ml 



]^' 21 '[lo, u>'] — 2ai[uj,uj']r[uj] l r[uj'] 1 + 9i>iai[uj,uj'](3i[uj,uj']r 2 [Lo]r 2 {io']r[uj] 1 l r{u>'} [ 1 
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l>0 

(21+1,21+1) r /i , n r n r iZ r nl 

H21 [a;, a; J = — ai[Ui,ui \a2[u>,0J \r[uj\ r[uj \ 

9i>i = for I = and 6i>\ = 1 for I > 1. As these expressions are symmetrical in the exchange of the two parameters 
u) and u/ (this was not the case for the general expressions 21 of 7^), equations (24) lead to 

L 

(L,L), 1 ifr iLujq /""do;' ,. (L,L) r 1 /n 



« L J 2 7 2tt v 

The various terms Xff correspond to a motional force evaluated at a delay time close to a multiple of the time 

of flight q between the two mirrors. For example, the term xTj ^ describes the response of the force F\ to the motion 
qi evaluated at times much shorter than the roundtrip time. As it could be expected, it does not depend upon the 
presence of the second mirror and is the same as if the mirror 1 were alone (compare /4i'°^ k = 2ai[w,a/] with 
equation 22). The other contributions to \n correspond to the modification of the damping force due to the presence 
of the second mirror. They depend upon the reflectivities of the two mirrors and appear at time delays close to a 
multiple of the roundtrip time 2q. 

The terms X2i' L ^ describe the force exerted upon one mirror when the other one moves. They depend upon the 
reflectivity of both mirrors and appear at time delays close to an odd multiple of the time of flight q. The lowest 
order term has a simple form 

jit^i' 1 ^ k u>'] = — oti[uj 1 uj']ot2[uJ 1 uj'] 

The moving mirror modifies the stress tensor of the vacuum field (modification described by the function a\). Then the 
radiation pressure experienced by the other mirror registers the modification of the stress tensor (detection efficiency 
described by the function a 2 ). 

In the limiting case of perfect reflection at all frequencies lower than w, one recovers the terms proportional to 
the third time derivatives in equations (23). The susceptibility functions diverge in this case but they are regular as 
soon as the mirrors have a small partial transmission. Indeed, the contributions corresponding to different numbers 
of roundtrips can be summed up to give in the resonance approximation 

R 2fl 1 [ W , W ']+a 1 [u), lI ;']fth W ']r 2 M r2 [a)'] e 2, ( u + u ')« 



7nM 
721 



a 1 [u},u)']a 2 [u,u']e i ( u+u ") q 



D[u,lu'} 

n 2»(w+w')9 



D[u,uj'] = 1 - r[w]r[w']e 

The denominator D characterizes the resonances of the Fabry Perot cavity considered as a mechanical resonator. 

Considering that the reflectivity coefficients may be approximated as constant functions from to w, the suscepti- 
bility functions are given by the simple expressions (r = rir 2 ) 

R n 4r i 

TnKwJ = 1 _ r2e2i(ljW ) g 

4 r 2 e i(w+a/)(} 
7*1 [W.W 7 ] = ~ t _ r 2 e 2»(u,W) g 

i ihLJ 3 r 2 e iuq 

X21 k 



6n 1 - r 2 e 2i "« 



corresponding to a motional force 

<5fl(f)> = ^ (rj (5q'{'(t) + r 2 Sq'{'(t - 2q) + rHq'(> ' (t - Aq) + . . ) - r 2 5q' 2 "(t - q) - rHq'^t - 3q) - . . .) 
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When compared with equation (23), one notes that the first time derivatives do not appear here because of the 
resonance approximation (the expression is valid only at high enough frequencies u ^> But the effect of imperfect 

reflection of the mirrors is now taken into account (r 2 and r\ are the reflection coefficients of the two mirrors for 
energy densities) and the susceptibilities are regular functions of the frequency. 



CONCLUSION 



The motional Casimir force constitutes a new type of interaction between two mirrors. As the stationary Casimir 
effect, it is associated with a modification of the vacuum stress tensor due to the field scattering upon the mirrors. 
But it is resonantly enhanced at the resonance frequencies of the optical cavity — . When compared with the case of 

r 2 

a single mirror, the enhancement may reach the value ._\ ^ . Consequently, the motional force might be very large 
p0[ with the high finesse cavities such as those used in cavity QED Q . 
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APPENDIX A: ANALYTIC PROPERTIES OF THE CORRELATION FUNCTIONS 

For the sake of clarity, we recall here the expressions of the susceptibility functions (see equations 2 and 19) 

Xi j (t)=2ie(t)^ j (t) (25a) 
€«(*) = Jt£J ^e-^-^M (25b) 
= (e(w) + e(w'))7«ku/] (25c) 

Since the susceptibilities are related to the retarded part of the correlation functions, their derivation relies upon the 
analytic properties of the coefficients 7^ . We show here how these properties can be inferred from the expressions of 
the coefficients 7^ in terms of the scattering coefficients, which are themselves analytic functions of the frequency (see 
equations 10). The coefficients 7y are obtained from products of two matrices Ti (see equation 16) which are functions 
of the scattering and resonance matrices S and R (see equation lib). Developping the corresponding expressions, one 
obtains rather lengthy expressions. 

However, these expressions may be simplified by using the following properties. First the S— matrix is unitary 
(S'^JS'lw]^ = /) and the matrix _R[w]i?[w]^ may be written in terms of the retarded (analytic for Imu > 0) and 
advanced (analytic for Imw < 0) components Q\u>\ and Qk^ (see equations 13). Then, it is also possible to reduce 
products _R[u->]<S'[lj]^ and S'^iJfo;]^ by noting that they determine the expression of the intracavity fields in terms of 
the output ones 

$ cav [H - R[ui[S[iS$ $outM 

Hence, i2[cj]5[w]T is an advanced response function and its adjoint S^wji?^]^ is a retarded one. Simple manipulations 
lead to 



S[w]R[u]i = R[uj] R[w]S[w]* = i?H f (26a) 

_L ( 52 

d[u] V sir 2 e iu "i si 



Using these properties, the coefficients 7^ are written as sums of terms, each of them being easily recognized as 
either a retarded or an advanced function of the two frequency parameters ui and to' 

7y k w'] = e iEj (Tr [P Ei P Ej + P- Si P- Sj } 

-Tr [P ei R[cu}P Ej R[u;'} T + P Ei R[^ P ej R{u'}*] 

- Tr [P^Rl^P-sjVf + P-eM^P-sjW] 

+ Tr [P £i (I + Q[u] + QH f ) P e] (I + Q[lo'] t + Q[u>']*)] 
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+ Tr [P_ £j (/ + Q[u] + QM 1 ) P- ej (I + Q[to'] T + Q[u/]*)] 
+e j( " W)9 Tr [P^S[iu}P £] S[u'} T ] + e - l(uJ+ul ' )q Tr [P^S^P^.^u/]*] 
-e' (uW) 'Tr [P_ £i P[w]P £ .PV] T + P- £l i?HP 3 i?K] T ] 
_ e -i(a>W) 9Tr ^.r^ p_ e _R[jy + p e .p[cj]tp_ e .S[o;']*] 

+e i(u,W)9 Tr [ P _ £i (J + Q[ w ] + Q[ w ]t) p £j (/ + Q[ W ']T + Q[a/]*)] 
+ e- 4 ^+ w ')« Tr [P £j (7 + + Q[a;]t) P_ £ . (7 + Q[^'] T + Q[u>']*)] 



Now, the susceptibility functions Xij are obtained by retaining the retarded terms and dropping the advanced ones 
in £y (see equations 25). Consider first the contribution of a term which contains the factor e(w), which is still an 
analytic function of a/ but not of oj. Its integration over ui' (see equations 25) provides a contribution to Xij{t) which 
is either retarded (vanishing for t < 0) or advanced (vanishing for t > 0). Some terms do not depend upon uj' and 
their integration contributes for a half to retarded terms and for a half to advanced ones as can be verified explicitly 
on expressions (25). The terms containing the factor e(o/) are computed in the same manner by considering the 
integration over lu. 

One finally obtains Xij(t) as an integral over two frequencies (see equations 2) with 



iTiLJLo' 



Xij K w 



( £ (c) 7 «[ W , W ']+ £ (c')75K,H) 



where 7^ is the retarded part of 7^ considered as a function of its second frequency parameter, that is (assuming 
q>0) 

7 £[w, "1 = e iej Q Tr [P £i P E . + P- £i P- E ,] 

-Tr [P £i PMP £3 P[ W '] T ] - Tr [P_jHP_ £ ^[ W f] 



•Tr 



Tr 



P £i (7 + Q[w] + Q[w]t)P Bj (~ + Q[u/] 



P_ £l (I + QH + Q[o;]t) P_ £ . ( - + Q[</] 



+e l(wW) « Tr [P_ £i S'[w]P £3 S*[w'] T ] 

_ e i{«+o/) a Xr [p_ £ .p[ w ]P £ .p[ w ']T + p_ £ .p[a,]P £ .p[a;'] T ] 
+ e < < w + w '>« Tr [P_ £i (7 + QH + Q[ w ]t) P £j (7 + Q[uj'] t )] 
+ e- i ( w+w ')« Tr [P £j (7 + QM + Q[a;]t) P_ £3 Q[ W '] T f 



Using the expressions (8), (9), (13) and (26) of the matrices S, R, Q and R in terms of the scattering coefficients 
describing the two mirrors, algebraic manipulations lead to the expressions (21) of the coefficients 7^. 

As discussed previously, the coefficients 7^ contain terms corresponding to retarded contributions and which appear 
also in the coefficients 7^. The other terms correspond to advanced contributions and have been dropped. A 
comparison between the two types of terms shows that the coefficients 7y can be deduced in a simple manner from 
the coefficients 7^ 



(27) 



This relation between the correlation function and the susceptibility function may be considered as the expression of 
the fluctuation dissipation theorem || for the present problem. With the help of expressions (21), it provides the 
explicit expressions of the coefficients 7^ in terms of the scattering coefficients. 



APPENDIX B: CONNECTION WITH SQUEEZING 



In the one mirror problem, it has been possible to compute the damping force by considering that the field scattering 
is modified when the mirror moves. It has been shown that this approach is completely equivalent to the linear response 
technique 0. This can also be shown in the two mirrors problem studied in the present paper. 
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For a single mirror, the modified scattering matrix Si can be written in a first order expansion in the mirror's 
motion 5qi around the position qi as 

5j[o;,a; ] = 27r5(o; — ui')Si[u] + SSi[u>,u>'} 

SSi[cj,u'} = iJ5q[uj - u r \e- iriu « (Si[u]r] - r}S i [u']) e lvuJ ' q ' 

The modification 5S and SR of the matrices associated with the Fabry-Perot can be derived from the elementary 
matrices SSi associated with each mirror. One computes 

<5S[w,c/] = ^P_ £i £S i [w,o/] (P £i + P_ Ei P[u/]) + ^P_ £i S4u;]P_ Ei <JP[a;,a/] 

z i 

z = 3 — i 

Then, one obtains the modified mean force through equations analogous to (11). The susceptibility functions given 
by the linear response theory are recovered at the end of lengthy calculations. 

This discussion shows that the motional Casimir force is connected to the problem of squeezing fig]. As a matter 



of fact, the motional modification of the field scattering corresponds to a squeezing of the input fields 15|,|7[. This 
squeezing generation requires that energy and impulsion be exchanged between the field and the mirrors. The motional 
Casimir force can be interpreted as a mechanical consequence of this effect. 

As in the single mirror problem, there exists an effective Hamiltonian which describes the squeezing effect (mod- 
ification of the field in response to the mirrors' motion) as well as the Casimir forces (mechanical action upon the 
mirrors in response to a variation of the field stress tensor). This effective Hamiltonian is the secular part (component 
at zero frequency) of the coupling (18) 

SH[Q} = Jdt SH(t) = J^J ^^^•^-^W Tr [^k^'] $ i I1 H^M T ] 
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